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We report the experimental reconstruction of a nonclassicality quasiprobability for a single-photon added 
thermal state. This quantity has significant negativities, which is necessary and sufficient for the nonclassicality 
of the quantum state. Our method exhibits several advantages compared to the reconstruction of the P function, 
since the nonclassicality filters used in this case can regularize the quasiprobabilities as well as their statistical 
uncertainties. A priori assumptions about the quantum state are therefore not necessary. We also demonstrate 
that, in principle, our method is not limited by small quantum efficiencies. 
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I. INTRODUCTION. 

The relation between classical optics and quantum optics 
is a fundamental topic in modern physics. The notion of non- 
classicality has been introduced by Titulaer and Glauber as the 
impossibility to describe optical field correlations of a specific 
state of light in terms of classical electrodynamics [l ]. There- 
fore, nonclassical states are the prerequisite of quantum ef- 
fects, which are of great interest in quantum optics and quan- 
tum information. 

The definition of nonclassicality is based on the Glauber- 
Sudarshan phase-space representation of a quantum state p, 

J d^aP{a)\a){a\, (1) 

where \a) denote the well-known coherent states lElO. In 
general, P{a) is a quasiprobability. If it has the properties 
of a classical probability distribution, the state p is a statis- 
tical mixture of coherent states, which are closely related to 
the classical behavior of the oscillator Conversely, a state 
is referred to as nonclassical if the P function shows some 
negativities. However, for many states, already for the single- 
photon state, this quantity is highly singular. Only in a few 
cases, when P{a) does not show singularities, it may be ac- 
cessible from experimental data. Single-photon added thermal 
states (SPATS) belong to such a class of states |4|, so that their 
P function could be reconstructed for some parameters |5 1. 

Due to the singularities of the P function, several different 
nonclassicality criteria have been developed. Simple inequal- 
ities often set bounds for classical states, which are violated 
for certain nonclassical states; we only mention bounds on 
moments, such as (higher-order) squeezing |6J, and probabil- 
ity distributions |7). However, violation of these criteria is 
only sufficient, but not necessary for verifying nonclassicality. 
On the other hand, there are complete hierarchies of criteria, 
often based on matrices of moments |8 1 or characteristic func- 
tions f9', TOl. However, nobody is able to check such hierar- 
chies completely, and practical application becomes involved 
for large matrices. 

In 111 II . a novel approach for nonclassicality detection was 



developed, which was based on phase-space methods. It has 
been shown that for any nonclassical state, there exists a so- 
called nonclassicality quasiprobability which illustrates the 
nonclassical property by negativities. Moreover, it was shown 
that a family of nonclassicality distributions, parameterized by 
a real filter width w, enables one to decide whether a state is 
nonclassical or not. If a state is nonclassical, one can always 
find some finite filter width w, such that the nonclassicality is 
observable as a negativity of the corresponding nonclassical- 
ity quasiprobability. Our method is suitable for experimental 
application, since it incorporates correct handling of statisti- 
cal uncertainties. Moreover, it does not require precognition 
about the state. 

In the present paper we demonstrate, by experimental ap- 
plication, the capability of the method. We examine single- 
photon added thermal states, whose P function could be re- 
constructed for a sufficiently large mean thermal photon num- 
ber Q. Here we overcome the problems occurring for ar- 
bitrary mean photon numbers, and therefore demonstrate the 
universality of the method of nonclassicality quasiprobabili- 
ties. 

The paper is structures as follows: In section II, we briefly 
review the approximate reconstruction of the P function of 
a SPATS and its limitations for relatively small mean photon 
numbers. In section III, we discuss the concept of nonclas- 
sicality quasiprobabilities and apply it to experimental data. 
Eventually in section IV, we consider the role of the quantum 
efficiency on the detection of nonclassicality. A summary and 
some conclusions are given in section V. 



II. APPROXIMATE RECONSTRUCTION OF 
P FUNCTIONS 

Let us briefly recall the reconstruction of a Glauber- 
Sudarshan representation as presented in Q. The starting 
point of the discussion was the characteristic function $(/?) 
of a SPATS, 

$(/?)= [!-(! + n)|/3p]e-«l^l^ (2) 
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By Fourier transform the resulting P function is derived as 

P(a) = 4^[(l + n)|ap-f^]e-l"l'/^ (3) 

which is a regular function. 

Experimental data were used for a mean photon number of 
n — 1.11 and a quantum efficiency of rj = 0.60. The func- 
tion was readily obtained from measured quadratures 
via 

It was directly sampled from 10^ data points. We observe 
that the experimentally obtained curve tends to zero within a 
fraction of its standard deviation, 

{*(/?)} -^p|'-i<f(/3)n- (5) 

However, the latter grows exponentially with |/3p. In order to 
calculate the Glauber-Sudarshan P function via Fourier trans- 
form, 

P{a) ^^J <fpe^f>'-"'P^{P), (6) 

one has to regularize ^(/S). In our previous work, we simply 
cutoff for |/?| > \Pc\- This is equivalent to the mul- 
tiplication of $(/3) with a rectangular filter, ilrcct(/3)i with 
^rcctiP) = 1 for |/3| < \(3c\ and n.^ctiP) = elsewhere. 
However, this method can only be applied in special cases. 
First, the state must be described by a well-behaved P func- 
tion, and its characteristic function has to decay sufficiently 
fast in order to justify the cutoff regularization. If this is not 
the case, one cannot perform the Fourier transform to obtain 
a P function. Second, the systematic error has to be esti- 
mated by some assumptions on the behavior of the charac- 
teristic function for large /3. We used the theoretical expec- 
tation of the characteristic function with properly chosen pa- 
rameters. For a completely unknown state, such a procedure 



becomes meaningless. In the case of the SPATS, we obtained 
the P function for mean thermal photon numbers of n > 1, 
for details see |5 1. 

For smaller mean thermal photon numbers the reconstruc- 
tion of the P function faces severe limitations. The smaller 
n, the broader the characteristic function becomes, and the 
larger is the statistical uncertainty at a reasonable cutoff pa- 
rameter. Therefore, one cannot find a useful trade-off between 
the systematic and the statistical error The former increase 
with lower and the statistical uncertainty is growing with 
larger |/3c|. For an example we refer the reader to the end of 
the next section. Under such circumstances, other nonclassi- 
cality criteria can be applied ifTOl . which are sufficient but not 
necessary. 



III. NONCLASSICALITY QUASIPROBABILITIES 

From a more general perspective to be used in the follow- 
ing, we may multiply the characteristic function by a filter 
function r2(/3), 

<^n{P) = ^[miP), (7) 

for a general study of such a scenario cf. lfT2l . The special 
case discussed so far is contained in this approach for a rect- 
angular filter. To obtain a quasiprobability - including the full 
information on the quantum state under study - by Fourier 
transform of the filter must not have zeros anywhere 

in the complex plane. 

In view of the radial symmetry of our quantum states, the 
Fourier transform of $s-j(/3) is given by 

Po(«) = - / bU2h\a\)^n{b)dh. (8) 
This function, together with its variance. 



can be readily calculated from the set of data. This expression 
for the variance of the quasiprobability is readily derived from 
the statistical sampling of the characteristic function accord- 
ing to Eq. Q. 



I 

A. Concept of nonclassicality quasiprobability 

Now we make use of the recently introduced concept of 
nonclassicality quasiprobabilities ifTTI . We introduce a so- 
called nonclassicality filter, = fiuj(/?), with the follow- 
ing properties: 

1. decays faster than cxp(— |/3p/2) for any filter 

width It; > in order to regularize the P function and 
its statistical uncertainty for any quantum state. 
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2. The Fourier transform of ilw{(3) is non-negative, such 
that negativities in the nonclassicality quasiprobabil- 
ity are unambiguously caused by the negativity of the 
state's P function. 

3. The parameter w scales the filter ^w{f3) such that for 
w oo the filter approaches one. Practically, we im- 
plement this by f]i(0) = 1, f^^(/3) = ^li^/w). 

4. The support of ilw{l3) is the complex plane, such that 
Eq. (j?]) is invertible for every (3. 

The first requirement ensures that the integrals in Eq. (|8]l 
and Eq. (j9]) are finite, when the filter is identified with 
the nonclassicality filter r2u,(/3). The second condition makes 
sure that the negativity of Pn{o:) unambiguously represents 
the nonclassicality of the state. In contrast to this, a rectangu- 
lar filter Street has a nonnegative Fourier transform. Hence, 
for such a filter negativities of (a) must not be interpreted 
as nonclassical effects without assumptions about the influ- 
ence of the filter on the regularized P function. In ||5], this 
made the estimation of systematic errors necessary. The third 
condition can be used for maximizing the significance of the 
observed nonclassical effects. On the one hand, a larger fil- 
ter width w may increase the negativities in the regularized 
P function, one the other hand, this will definitely increase 
the variance. We may tune w in order to optimize the statisti- 
cal significance S of the negativities, defined via 



Pnja) 



(10) 



The fourth requirement is of fundamental importance: It en- 
sures that the regularized function Pq (a) still contains the full 
information on the quantum state. Such Pq functions, which 
fulfill all the four conditions, we refer to as nonclassicality 
quasiprobabilities. 

A proper nonclassicality filter can be constructed by the au- 
tocorrelation of a function 



(11) 



with the normalization U — j d^/3'|a;(/3')p. The width pa- 
rameter is introduced by ^yj{fi) = Vli{l3/w). The positivity 
of its Fourier transform is guaranteed by the properties of any 
autocorrelation function. Furthermore, if w(/3) decays suffi- 
ciently fast, as required by condition (1), does as well. 
Therefore, by choosing 



(12) 



and calculating the autocorrelation, we obtain a suitable rep- 
resentative of a nonclassicality filter 



B. Experimental preparation of SPATS 

The single-photon-added thermal states are generated in a 
conditional way by exploiting the parametric amplification at 



the single-photon level in a nonlinear type-I /3-barium borate 
(BBO) crystal pumped by the second harmonic of a mode- 
locked picosecond Ti:sapphire laser (see Fig.[T]l. 

When no extra field is injected in the crystal, a pump photon 
can be converted into two spontaneously and simultaneously 
generated photons (named signal and idler) correlated in fre- 
quency and in momentum. The click of the on/off avalanche 
photodetector D, which is placed in the idler path after narrow 
spectral and spatial filters (F), is used to conditionally prepare 
a single photon in a well-defined spatiotemporal mode of the 
signal channel 1.13. ,14J . On the other hand, if the parametric 




Figure 1. Scheme for the conditional excitation of a thermal light 
state (denoted by pin) by a single photon. A click in the on-off detec- 
tor D prepares the photon-added thermal state pout and triggers its 
balanced homodyne detection (BHD). 

crystal is seeded with some light, described by the operator 
/5i„, stimulated emission comes into play, and single-photon 
excitation of such a state, always converting it into a non- 
classical one, is conditionally obtained when one photon is 
detected in the idler mode ifTSl [161 . Field quadratures of the 
output signal state are then conditionally measured on a pulse- 
to-pulse basis using an ultrafast balanced homodyne detection 
scheme [17]. 

Here we used a pseudo thermal source, obtained by insert- 
ing a rotating ground glass disk in a portion of the laser beam 
(see ifTOl ITSi 19 1), for injecting the parametric amplifier and 
producing SPATS. 



C. Experimental nonclassicality quasiprobabilities 

To illustrate the power of nonclassicality quasiprobabilities, 
let us consider a SPATS with n = 0.49 and rj = 0.62. In 
Fig.|2] we show the experimentally reconstructed characteris- 
tic function $oxp(/3)- Obviously, the Fourier transform of this 
quantity does not exist as a regular function, since $oxp(/?) 
does not approach zero for large /? as its theoretical expecta- 
tion $th(/3) does. This is due to the fact that the uncertainty 
grows exponentially. Although both the experimental result 
and the theoretical expectation agree within two standard de- 
viations, the former may even diverge within the divergent 
noise level. In contrast to the results in ||5J, it is not obvious 
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Figure 2. Characteristic functions of a SPATS with n — 0.49 and 
ri = 0.62: the experimental result <l>oxp(/3), the theoretical expec- 
tation $th(/3), and the result ^niP) of filtering for the filter width 
w = 1.4. The shaded area corresponds to one standard deviation. 




Figure 3. Nonclassicality quasiprobability of a SPATS for the same 
parameters as in Fig. |2] The blue shaded area corresponds to one 
standard deviation. 



just from experimental data, that the characteristic function 
tends to zero for large /3. Moreover, it is not possible to find 
a reasonable compromise between a low systematic error (re- 
quiring a large cutoff parameter \f3c\), and a statistical uncer- 
tainty being sufficiently small to obtain significant negativities 
in the filtered P function. 

The application of the nonclassicality filter with a width 
w — lA leads to an integrable characteristic function <I>n(/3). 
We emphasize that the shown curve is obtained from the ex- 
perimental data. We also calculated its standard deviation, 
which is included in the line thickness. Therefore, this func- 
tion is suited for deriving the corresponding nonclassicality 
quasiprobability by Fourier transform. Figure [3] shows the re- 
sult. We observe a distinct negativity at the origin of phase 
space, with a significance of 15 standard deviations. By the 
definition of nonclassicality quasiprobabilities, this negativity 
is solely due to the nonclassicality of the state. 

We may also reconsider the data of the SPATS with n = 
1.11 by our filtering procedure. We obtain a nonclassical- 
ity quasiprobability PniP) which looks similar to the one in 
Fig. [3] By optimizing the filter width to w — 1.3, we get a 



maximum significance of 7.6 standard deviations, which ex- 
ceeds that for the rectangular filter. More importantly, the es- 
timation of a systematic error - which was previously needed 
for the assessment of nonclassicality and had been based on 
a-priori theoretical assumptions - now becomes superfluous. 

D. Limits of rectangular filtering 

For comparison, let us try to regularize the state with n = 
0.49 and rj ~ 0.62, whose characteristic function is shown 
in Fig. |2] with a rectangular filter. Since this is not a non- 
classicality filter of the type defined in Sec. Ill A, we have 
to consider the systematic error, which comes from the cutoff. 
From Fig.|2] two possible cutoff parameter may be reasonable: 
On the one hand, we may choose |/3c| = 2.2, since from this 
point the deviations of the theoretical and experimental char- 
acteristic functions <i>t,h(/3) and ^cxpW) grow strongly. The 
corresponding regularized P function with its standard devi- 
ation (shaded blue area) and systematic error (dark red area) 
is shown at the left side of Fig.|4] We clearly see that the sta- 
tistical uncertainty is negligible, while the systematic error is 
partly even larger than of the P function. Therefore, a-priori 
assumptions about the state, which are the basis for the es- 
timation of the systematic error, are crucial. Moreover, the 
systematic error is larger than the size of the negativity. 
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Figure 4. (color online) P functions regularized with a rectangular 
filter. The cutoff parameter at the left side of the figure is rather 
small, leading to large systematic errors (red shaded). For a larger 
cutoff parameter, the statistical uncertainty (blue shaded) becomes 
dominant, see the right side of the figure. 

On the other hand, we may set |/3c| = 3.8, where the exper- 
imental characteristic function $cxp(/3) is close to zero. The 
resulting P function is shown at the right side of Fig.|4] In this 
case, the systematic error is small, while the statistical uncer- 
tainty is dominating. The significance of the negativity of the 
filtered P function is less than 1.4 standard deviations. 

Therefore, we are not able to obtain statistically significant 
negativities by applying a rectangular filter, if we simultane- 
ously try to achieve low systematic eiTors due to the regular- 
ization. The reason is that the characteristic function of a state 
with such a small mean photon number approaches zero only 
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for larger /?, where the standard deviation is already extremely 
large. This problem can only be overcome by using suitable 
nonclassicality filter, which do not require the consideration 
of systematic errors in order to verify nonclassicaUty. 



IV. INFLUENCE OF THE QUANTUM EFFICIENCY 

Let us assume that the characteristic function $(/3) de- 
scribes a quantum state, but we measure this stale with a non- 
unit quantum efficiency r]. The characteristic function of the 
measured state is given by 

$(^;r?)=$(V^/3). (13) 

Since the corresponding P function is the Fourier transform, 
this rescaling of /? leads to a rescaling of the argiunent of the 
P function, 

''^"^''^ = ^'^(^)- '''^ 

Therefore, a non-unit quantum efficiency does not degrade 
nonclassicality, and removal of losses from experimental data 
only rescales the P function, but does not affect its course 
and uncertainty. In contrast, this is not the case for the Wigner 
function whose negativities more and more disappear with de- 
creasing ry-values. 

How is the situation for nonclassicality quasiprobabilities? 
For answering this question, we note that the filtered charac- 
teristic function of a state, which suffered losses, reads as 

$n(/3;r?) = = (15) 

We observe that the characteristic function of the lossy state, 
^n{P', v)' ™d the one of the ideal state, ^n{0: 1) are not con- 
nected by simple rescaling. However, they can be easily con- 
nected if one allows to rescale the filter width as well: 

Hence, the characteristic function of a lossy state with some 
width w, which is on the left side of the equation, is given 
by the characteristic function of the ideal state with filter 
width by additionally rescaUng the argument /3 
Removing the losses by post-processing and simultaneously 
adapting the width of the nonclassicality filter leads simply to 
scaled results with no better significance. Therefore, the re- 
moval of losses does not uncover nonclassical effects which 



do not already appear in the nonclassicality quasiprobabilities 
of the lossy state. 

Finally, we note that our method may visualize nonclas- 
sicality even for rather small quantum efficiencies. Simula- 
tions show that we find for the SPATS with n = 0.49, for our 
sample of 10^ data points, negativities in the nonclassicality 
quasiprobabiUty with a significance of at least three standard 
deviations, if > 0.4. This ?7-value is only limited from be- 
low by the size of the sample. In contrast, for 77 < 0.5 the 
Wigner function is always nonnegative. Therefore, the nega- 
tivities of the nonclassicality quasiprobabilities are more sen- 
sitive to nonclassical effects than negativities of the Wigner 
function. 

V. SUMMARY AND CONCLUSIONS 

We have applied the concept of nonclassicality quasiproba- 
bilities to experimental data of single-photon added thermal 
states. Even though the Glauber-Sudarshan P function of 
these states is regular in general, its approximate reconstruc- 
tion is feasible only for a certain parameter range. Moreover, 
it requires to make use of some precognition on the state under 
study. Our quasiprobability approach does not require such a 
precognition, it works for any quantum state - even when the 
P function is strongly singular - and it also suppresses the 
experimental sampling noise. 

We have shown that the nonclassicality filters needed in our 
procedure, suppress the exponential growth of experimentally 
determined characteristic functions, which yields integrable 
functions. Hence, Fourier transform delivers nonclassicality 
quasiprobabilities with finite statistical uncertainties. By op- 
timization of the filter width, significant negativities in the 
quasiprobabilities are found for nonclassical states. For our 
chosen example, an approximate reconstruction of the P func- 
tion was shown to be no longer useful. 

With the accessible set of 10^ data points we could demon- 
strate negativities in the experimentally determined nonclas- 
sicality quasiprobability with a significance of about 15 stan- 
dard deviations. This result is solely limited by the statistical 
uncertainties caused by the finite size of the available set of ex- 
perimental data, it could be further improved by extending this 
set. We have also considered the role of imperfect detection. 
In fact, the detection efficiency can be completely eliminated 
by a proper rescaling of arguments in our functions. Thus, 
even with a rather small efficiency one can identify all non- 
classical effects, provided the sampling noise is sufficiently 
small. 
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